Temperature dependence of the contact in a unitary Fermi gas 
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The contact T, introduced by Tan, has emerged as a key parameter characterizing universal 
properties of strongly interacting Fermi gases. For ultracold Fermi gases near a Feshbach resonance, 
the contact depends upon two quantities: the interaction parameter l/ikpo), where kp is the Fermi 
wave-vector and a is the s-wave scattering length, and the temperature T/Tp, where Tp is the Fermi 
temperature. We present the first measurements of the temperature dependence of the contact in a 
unitary Fermi gas using Bragg spectroscopy. The contact is seen to follow the predicted decay with 
temperature and shows how pair-correlations at high momentum persist well above the superfiuid 
transition temperature. 
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Two component ultracold Fermi gases near Feshbach 
resonances provide an archetypal setting to explore uni- 
versal behaviors [TH3]. Universal systems should satisfy 
two requirements: firstly, the gas must be dilute enough 
that the mean interparticle spacing rf^/^ greatly exceeds 
the range of the interaction potential ro, and, secondly, 
the interactions, characterized by the s-wave scattering 
length, a, should be sufficiently strong that a greatly ex- 
ceeds n~^/^. All Fermi systems that satisfy these require- 
ments will behave identically on a scale given by the aver- 
age particle separation, independent of the details of the 
interaction potential. Recent theoretical work by Tan Pl- 

and others [TUTT] has identified several exact relations 
applicable to Fermi systems in the universal regime. The 
central parameter in these relations is a quantity called 
the contact I, which forms a link between microscopic 
and macroscopic system properties. 

Contact quantifies the likelihood of finding two inter- 
acting fermions with very small separation and is closely 
linked to the pair-correlation function [4 . In strong anal- 
ogy with the phase diagram of the Bose-Einstein con- 
densate (BEC) to Bardeen-Cooper-Schrieffer (BCS) su- 
perfiuid crossover of two-component Fermi gases [12 , T 
depends upon two parameters: the dimensionless inter- 
action strength l/{kpa), where kp is the Fermi wave- 
vector, and the relative temperature T/Tp, where Tp is 
the Fermi temperature. Previous theoretical |7l [13] and 
experimental [Tl'-TOI work has investigated the interac- 
tion dependence of the contact, and a number of recent 
studies have calculated the temperature dependence of 
the contact [TTHT^ . To date however, there have 
been no measurements of this latter dependence. 

In this letter, we report the first measurements of the 
temperature dependence of the contact using Bragg spec- 
troscopy of a ^Li Fermi gas at unitarity. Bragg spec- 
troscopy allows for quantitative measurements of the 
static structure factor S{k) which is directly proportional 
to the contact at high momenta. We extract the first 



and second moments from our Bragg spectra and use 
these to obtain the dynamic structure factor S{k, lo) and 
from this S{k) and the contact. Our results are in good 
agreement with theoretical predictions and indicate that 
pair-correlations at high momenta persist well above the 
critical temperature for superfiuidity. 

Tan's exact relations for Fermi gases near the BEC- 
BCS crossover marked a dramatic development in the 
understanding of highly correlated Fermi systems [3]- 
|6j. Strongly interacting Fermi gases represent a diffi- 
cult theoretical challenge as the large scattering length 
leaves no small parameter that can be used in pertur- 
bative theories. Having invoked the contact. Tan was 
able to derive a number of relations linking the micro- 
scopic parameters to macroscopic properties such as the 
energy and momentum distribution. Contact is defined 
as I = lim/j_>.oo fc*«-o-(fc) where the momentum density 
of a particular spin component n„{k) decays with 1/k^ 
at large k. Thus I quantifies the amplitude of this high 
momentum tail. A more intuitive understanding of I is 
evident through its relation to the two-body correlation 
function between spin-up (f) and spin-down {],) fermions 

valid for tq < r < kp^ , where the contact appears as a 
prefactor |3]. The contact therefore quantifies the likeli- 
hood of finding two fermions at distances small compared 
to the many-body length scales. 

Contact is closely linked to the pairing temperature 
T* in the high momentum limit. A number of experi- 
ments have measured the contact using photo-association 
[3 [2] , radio-frequency spectrocopy [15j , the tail of the 
momentum distribution [15] and Bragg spectroscopy [TB] . 
These measurements demonstrated the predicted decay 
of the contact through the transition from the BEC to 
BCS sides of a Feshbach resonance. Pairing (T*) in the 
phase diagram of the BEC-BCS crossover is also strongly 
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temperature dependent [T2]. Our measurements demon- 
strate the build up of pair-correlations from T ^ T* 
down to T — > 0; in the momentum range where Tan's 
relations hold. 

The Fourier transform of gj^ (r) yields the static struc- 
ture factor S^i{k) EOJIIT], which can be quantitatively 
measured using Bragg spectroscopy [IS]. In the limit 
where a ^ oo the second term in Eq. 1 vanishes and the 
S^l{k) is given by 

(2) 

where I/{Nkp) is the dimensionless contact and N is 
the total number of atoms. At large momenta, k ^ kp, 
the total static structure factor S{k) = 1 -|- S^i{k) as 
the spin-parallel component approaches the uncorrelated 
value of unity at large k PTI 1^ . Therefore we can read- 
ily determine the spin-antiparallel component. by 
measuring the total structure factor. 

Experimentally we create a unitary Fermi gas by evap- 
oratively cooling a balanced mixture of ^Li atoms in the 
two lowest spin states \F = 1/2, mp = ± 1/2 ) at a 
magnetic field of 834 G in an optical dipole trap. After 
transferring the atoms into a second deep dipole trap, 
we obtain N^r = 170, 000 atoms in each spin state at a 
temperature of 0.09 ± 0.03 T/Tp. The trap frequen- 
cies in the final trap are {u!x,ujy,0Jz)/2TT = (24.5,65,230) 
Hz, (w = 71.5 Hz) giving a Fermi energy Epl{2nh) = 
7.2 kHz, where Ep = huj{ZNf^^ and kp = ^/2mEp/h 
= 2.9 fim^^ . At 834 G the s-wave scattering length di- 
verges and collisions are unitarity limited. 

We vary the temperature of the cloud by suddenly 
switching off the dipole trap for a variable time before 
quickly ramping it on and holding the cloud for 400 ms 
to rethermalize In this way we can re- 

peatably heat the cloud to temperatures up to 1.1 T/Tp 
without loss of atoms. To determine the temperature, 
we image clouds at unitarity after 2 ms expansion and 
obtain an empirical temperature T by fitting Thomas- 
Fermi profiles to the images |23J. While the conversion 
from T to the true temperature T/Tp = Ta/1 -|- /3, where 
(3 K, —0.58 is the universal parameter [21], is not exact, 
we have independently calibrated it through measure- 
ments of the mean energy per particle and compared 
these to predictions based on a Noziere-Schmitt-Rink 
(NSR) theory |i25) which agrees well with thermodynamic 
measurements of the equation of state |26If28| . The two 
methods agree to better than 0.05 Tp over the full range 
of temperatures studied here. 

Bragg scattering is achieved by illuminating a cloud 
of atoms with two laser beams that have a small fre- 
quency difference w and intersect at an angle 9. This 
creates a periodic potential which moves at a velocity 
uj/k, where k = 47r sin (0/2)/A and A is the wavelength 
of the Bragg light. Measuring the response of the cloud 



to a sequence of Bragg pulses as a function of w, yields a 
Bragg spectrum from which quantitative information on 
the dynamic structure factor can be obtained |29| . 

To perform Bragg spectroscopy, we use two laser beams 
that intersect at 6* = 49.5°. At this angle, k = 2.7kp and 
the resonant recoil frequency for Bragg scattering of free 
atoms is Wr/(27r) = 51.6 kHz. The Bragg lasers are de- 
tuned ~2.5 GHz from the scattering transition to avoid 
spontaneous scattering and the Bragg pulse duration is 
200 /iS. For our beam intensities this duration is well be- 
low the two-photon Rabi cycling period ensuring spectra 
are obtained in the linear response regime. Each of the 
two spin states in the ^Li gas couple almost equally to 
the Bragg lasers to within 4%. We measure the resulting 
atomic distribution, n{x, y), a further 3 ms after applying 
the Bragg pulse. The trap laser is switched off immedi- 
ately after the Bragg pulse. Both Bragg scattering and 
imaging take place at 834 G. From these images we ob- 
tain line profiles n(x) by integrating over the y-direction 
perpendicular to the Bragg scattering. 

We perform a series of experiments to acquire a se- 
quence of profiles n{x) as the frequency difference be- 
tween the two Bragg lasers w is varied. Quantitative 
analysis is achieved by evaluating the first and second 
moments (x) and {x'^), respectively, for every line pro- 
file, where (a:™) = J2i ^T''^i^i)/ J2i "-(^^i) ^-^d the sum is 
over all pixels i. At large w no Bragg scattering occurs 
and the moments (xq) and (xg) provide a reference center 
of mass position and mean square cloud width, respec- 
tively. The difference between these reference moments 
and those obtained when Bragg scattering occurs pro- 
vides two ways to quantify the effect of the Bragg pulse. 
The center of mass displacement due to the Bragg pulse 
is given by AX(uj) = {x{uj)) — (xq), which is proportional 
to the momentum transferred to the cloud, and the in- 
crease in mean square cloud size Act^(w) — {x^{uj))~{xl) 
is proportional to the energy transferred Spectra ob- 
tained the the first and second moments, at temperatures 
of T/Tp = 0.10, 0.29, 0.58 and 1.0, are presented in Fig. 
1(a) and (b), respectively. Each spectrum is an average 
of four individual spectra. 

As the duration of the Bragg pulse is relatively long, 
the effects of Fourier broadening on our spectra are small 
and we ignore them in our analysis. For a given run of the 
experiment, with a particular value of w, the momentum 
{energy} transferred to the cloud will be proportional 
to k ■ [S{k,uj) - S{-k, -w)] {uj ■ [Sik,Lo) - S{-k, -uj)]} 
|29| . Apart from the factors of k and to, the pre- factors 
multiplying the difference of the positive and negative 
components of the dynamic structure factor are iden- 
tical, so the energy transferred is equal to w/k times 
the momentum. We include the negative component 
S{—k, — w) in our analysis as this describes de-excitation 
from high lying states which can become significant at 
the higher temperatures and momentum (fc = 2.1kp) we 
consider. To account for this, the principle of detailed 
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Figure 1. Bragg spectra and dynamic structure factor of a 
unitary Fermi gas. Measured change in (a) the first moment 
(center of mass displacement) AX (to) and (b) the second mo- 
ment (increase in mean square size) Aa^{uj) as a function of 
the Bragg frequency uj at temperatures T/Tf = 0.10, 0.29, 
0.58 and 1.0. (c) Dynamic structure factor S{k,uu) obtained 
from a weighted average of the first and second moments. 
Points are the experimental data and the solid lines are dou- 
ble Gaussian fits normalized according to the /-sum rule. 

balancing is employed which states that S{—k, —lo) — 
exp {—-i^Hp)S{k,oj) where fc^ is Boltzmann's constant. 
Thus both the first and second moments provide a mea- 
sure of S{k,ijj), as k,u} and T are all known. 

In Fig. 1(c) we show the experimentally determined 
dynamic structure factors obtained from a weighted av- 
erage of Ax{uj)/k and Aa'^{uj)/uj, for the same tem- 
peratures. The weighting of each contribution is in- 
versely proportional to the relative uncertainty in the 
respective measurement points. The detailed balance 
term [1 — exp (— ^^^)] has also been divided out to yield 
S{k,uj). This increases the noise in the data at low fre- 
quencies where this term becomes small, particularly at 
higher temperatures. Along with the data points are 
fits based on two Gaussian functions, centered near uJr/2 
and LOr to account for pair and free atom excitations in 
S{k,uj), respectively. While the true structure factors 
will not necessarily be well described by two Gaussians 
|30j . we have estimated that the errors introduced in ap- 
proximating S{k,ijj) with two Gaussians is at the level 
of a few percent for the coldest clouds and less at higher 
temperatures. The spectra are normalized according to 
the /-sum rule |3T], such that the fitted functions for 
S{k,uj) satisfy the integral h J LdS{k,uj)duj = Ncor- At 
low temperatures S{k,uj) is dominated by excitations at 



pair frequencies uJr/2 and the peak in S{k,uj) shifts to- 
wards ujr as T approaches Tp. The change from the low 
T to high T limits occurs smoothly over the temperature 
range covered as expected from quantum virial expansion 
calculations At high temperatures S{k,uj) contains 
significant weight at negative frequencies. 

From these dynamic structure factors, we can obtain 
the static structure factor S{k), defined by NS{k) = 
h J S{k, ui)duj. As we do not measure S{k, uj) at negative 
frequencies, and because of the enhanced noise in the ex- 
perimental data at low frequencies, we use the integral 
of the double Gaussian fits to S{k,u!) over uj to obtain 
S{k). As these satisfy the /-sum rule the integral pro- 
vides a robust measure of S'(fc) [T2]. Equation (2) allows 
us to directly link S^i{k) = S{k) — 1 to the contact. Ex- 
tracting the contact from the dynamic structure factors 
(Fig. 1(c)) for all measured temperatures gives the data 
shown in Fig. 2. The vertical error bars are the statistical 
uncertainty based on the range of values obtained from 
the individual spectra and the horizontal error bars in- 
dicate the uncertainty in the temperature measurement. 
Also shown in this figure are different theoretical calcula- 
tions for I. The three full lines are obtained from strong- 
coupling theories based on the many-body t-matrix ap- 
proximation: a non-self-consistent (GqGo) theory (pur- 
ple) [13], a self-consistent (GG) theory (green) [18j and 
a Gaussian pair fluctuation theory (GPF) (brown) |19| . 
The dashed gray lines show the results of a quantum 
virial expansion to second (long dashed) and third (short 
dashed) order |19j . 

Calculations of the zero temperature contact typically 
predict values between 3 and 3.4 |13| [TBI UHl HI] con- 
sistent with our measurement at the lowest temperature 
3.11 ± 0.23 and an earlier value of 3.40 ± 0.18 extracted 
from measurements of the equation of state [33]. At 
higher temperatures, the contact is seen to decay mono- 
tonically from ~ 3 down to below 0.5 over the range 
0.1 — ITp, in good general agreement with theory. The 
different theoretical methods each predict slightly differ- 
ent behavior near the critical temperature for superfluid- 
ity Tc « 0.2 Tp. The uncertainties in our measurements 
and limited number of data points prevent us from iden- 
tifying one theory as more accurate than another. With 
further improvements, and a closer study of the region 
near T^, it should be possible to resolve the discrepancies 
between different theories. 

We note that X is signiflcant at temperatures well 
above Tc and shows the gradual build up of pair- 
correlations below Tp. At the momentum we have stud- 
ied {2.7hkp) the 1/fc^ momentum tail dominates the pair- 
correlation function so these measurements should not 
be taken as evidence for pseudogap pairing. Whether a 
connection between short-range (r <^ kp^) and medium- 
range (r ^ kp^) pair-correlations can be found remains 
an open, but important, question [13] that may offer new 
insights into pairing both below ^34il35| and above Tc |36J. 
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Figure 2. Temperature dependence of the contact over the 
range T/Tp ~ 0.1 — 1.0. Data points are the measured values 
of I obtained from the static structure factors. Sohd lines are 
predictions based on three different strong-coupling theories 
as described in the text |13l 1181 119| . The long-dashed and 
short-dashed lines are second and third order quantum virial 
expansion calculations, respectively |19| . 



In summary, we have presented the first measurements 
of the temperature dependence of Tan's contact in a 
strongly interacting Fermi gas. These were achieved us- 
ing Bragg spectroscopy to make quantitative measure- 
ments of the dynamic and static structure factors. Our 
results indicate that the contact, and hence the high mo- 
mentum component of the pair-correlation function, re- 
mains significant at temperatures well above Tc, in good 
agreement with theoretical predictions. More extensive 
studies of the region near are needed to resolve the dis- 
crepancies between different theoretical approaches. Ex- 
tracting the homogeneous contact from experiments on 
trapped systems would allow a more sensitive probe of 
these discrepancies |13} IT9| . Bragg experiments at lower 
momentum (fc ~ kp) may also elucidate a connection 
between the contact and pseudogap pairing. 
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